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(57) ABSTRACT 


Method and system (i) to determine when a selected target 
molecule is present or absent in a fluid, (2) to estimate con- 
centration of the target molecule in the fluid and (3) estimate 
possible presence of a second (different) target molecule in 
the fluid, by analyzing differences in resonant frequencies of 
vibration of a thin beam suspended in the fluid, after the fluid 
has moved across the beam. 
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NANOSTRUCTURE SENSOR OF PRESENCE 
AND CONCENTRATION OF A TARGET 
MOLECULE 

This application claims priority based upon, and hereby 5 
incorporates by reference, a Provisional Application, No. 
60/849,494, filed 28 Sep. 2006. 

ORIGIN OF THE INVENTION 

to 

This invention was made, in part, by one or more employ- 
ees of the U.S. government. The U.S. government has the 
right to make, use and/or sell the invention described herein 
without payment of compensation, including but not limited 
to payment of royalties and other licensing fees. 15 

FIELD OF THE INVENTION 

This invention is related to sensing of presence of a target 
molecule in a fluid, using nanostructures. 20 

BACKGROUND OF THE INVENTION 

Several methods are presently being developed to sense 
presence of a target molecule having a small or modest con- 25 
centration in a fluid. However, estimation of the concentration 
of the target molecule is unavailable, where the fluid is sub- 
stantially motionless relative to the sensor, in part because the 
fluid in equilibrium will have substantially the same concen- 
tration everywhere. At least two questions should be 30 
answered here: (1) Is the target molecule present in the fluid; 
and (2) If the target molecule is present, what is its estimated 
initial concentration? 

What is needed is a nanotechnology approach that works 
with fluids (liquids and gases) to determine whether a target 35 
molecule is present in the fluid, and when the target molecule 
is present, to estimate an initial concentration of this sub- 
stance. Preferably, the approach should allow a distinction to 
be made, in appropriate circumstances, between presence of 
first and second different target mo lecules and should indicate 40 
when the target molecule is attached to, and has saturated, a 
sensor. Preferably, the approach should distinguish between 
presence of a given target molecule and presence of a mol- 
ecule that competes with the target molecule for attachment. 

45 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
nanostructure, having known mechanical properties and hav- 
ing one, two or more spaced apart beams partly coated with a 50 
selected chemical substance CS that binds to or has an affinity 
for the target molecule TM. The fluid is encouraged to move 
relative to the nanostructure and to interact with the chemical 
substance CS. In the absence of the target molecule TM, the 
microchannel has one or more predictable resonant frequen- 55 
cies of longitudinal vibration. Where the target molecule 
becomes attached to the beam the resonant frequency of 
vibration changes to another predictable value, according to 
one or more dispersion relations. Where two or more spaced 
apart beams are provided, the variation in number of target 60 
molecules TM attached to each of the two or more beams 
allows estimation of the concentration of the target molecule 
in the fluid. The invention allows the fluid to move relative to 
different sensor components to allow sensing of the presence 
of, and estimation of the concentration of, the target mol- 65 
ecule. From a knowledge of the maximum concentration of 
the target molecule TM attached to the beam(s), the invention 


2 

can also indicate whether the target molecule is saturated so 
that an estimate of the concentration in the fluid may be 
inaccurate. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates a channel constructed 
according to an embodiment of the invention. 

FIG. 2 is a flow chart of a procedure for practicing the 
invention. 

FIG. 3 is a flow chart of a procedure for practicing the 
invention. 

FIG. 4 illustrates a mechanism for determination of a wave- 
number or frequency solution according to the invention. 

FIG. 5 schematically illustrates reflection or diffraction of 
an incident electromagnetic wave by a vibrating beam. 

FIGS. 6A 6B and 6C illustrate systems for estimating 
concentration of a target molecule. 

FIG. 7 graphically illustrates variation of a parameter a 
with target molecule concentration r \ . 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

FIGS. 1 and 2 schematically illustrate an embodiment of 
the invention, using a sequences of one or more substantially 
identical, thin rectangular parallelepiped beams ll-n(n=l, . . 

. , N; Nis 1 ) 5 each attached at both ends to a non-movable wall, 
of solid material of dimensions 2L X by 2L y by 2L Z , with 
known material properties and having a longitudinal axis LL 
in the x-direction. The beam 11a selected chemical substance 
CS of width preferably 2L y applied in an orientation that is 
substantially along the longitudinal axis LL. A fluid (liquid or 
gas) 15 having a known density is allowed to flow across the 
beam 11 and the CS. Initially, the fluid 15 is absent, and one 
or more resonant frequencies of each reference beam before 
the fluid is introduced, is estimated. 

In one embodiment the chemical substance CS is a thin 
coating of a selected antibody, applied to vertically oriented 
nanostructures on a substrate, or applied directly to the sub- 
strate. The substrate plus CS coating provides a reference 
beam. Where two or more beams are provided, the beams are 
spaced apart in a selected fluid flow direction. The fluid 15 
may or may not contain an antigen or other target molecule 
TM that binds to the antibody, and presence or absence of the 
TM in the fluid 15 is to be determined. The chemical sub- 
stance CS need not be the full antibody: provision of a vari- 
able segment “tip” of one of the three arms of an antibody will 
often suffice to bind to its corresponding antigen, as discussed 
by W. T. Gibbs in “Nanobodies, ” Scientific American, August 
2005, pp. 79-83. Use of a variable segment tip of an antibody 
as the chemical substance CS, rather than the bulky antibody 
itself, allows distribution of a larger number of these tips on 
the nanostructure surfaces, and thus allows more robust 
attachment to the antigen. Extremes of pH and/or temperature 
are less likely to denature the tips, as compared with denatur- 
ing of the lull antibody. 

In the absence of the target molecule TM, each of the 
beams 11-n is homogeneous, has a known mass w per unit 
length, and satisfies a relation 

(d/dx)(EA du/dx)+(d/dy)(EA du/dy) + (d/dz) (EA du/dz)= 

(wA)(d 2 u/dt 2 ), (1) 

for longitudinal beam vibration, where E and A are Young’s 
modulus and cross sectional area associated with the beam, 
respectively, and u(x,y,z,t) measures vibration distance in the 
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x-direction of a beam volume containing the location (x,y,z). 


Suitable boundary conditions are 

u{x=±L x , y,z,t)= 0, (2) 

The unconstrained solutions for this defining relation for the 

u{x,y,z,t )= sin {A^jsin {k y y] sin {kyz) cos {( 0 0 (f-f 0 )}, (3) 

(D 0 (ref ) 2 =(E/w){k 2 +k y 2 +k z 2 }, (4) 

k=2xnJ2L x , (5-1) 

k=2 zmJ2L y , (5-2) 

k=2itnJ2L z , (5-3) 

n x ,Hy, ^=1,2,3,... (6) 


The unmodified, homogeneous beam (coated with CS over a 
limited region, -L 1 ^x^L 2 ) will vibrate longitudinally (x-di- 
rection) with a sequence of resonant angular frequencies 
w^n^n^nj given by Eq. (4). The added mass for CS is 
usually negligible and is ignored here. 

Where one of the target molecules TM becomes attached to 
the beam 11 in a second region (-L 1 ^x^L 2 ), modified solu- 
tions for a first region (-L x ^x^-L 1 ) and for a third region 
(L 2 ^x^L x ) are formally the same as in Eqs. (3) and (4), but 
the wavenumber parameters k x , k y , and k z , need not be the 
same for the first, second and third regions and are not nec- 
essarily determined by Eqs. (5-1), (5-2) and/or (5-3). With the 30 
target molecule attached to the second region, the defining 
equations become 

d 2 u 1 '!dx 2 +d 2 u 1 'Idy^+d^tl 1 Vdz 2 ^ w/E) (< d 2 u 1 '/dt 2 ), (7-1) 

d 2 u2 , /dx 2 +d 2 u2'/dy 2 +d 2 u2'/dz 2 =(w , /E)(d 2 u2'/dt 2 ), (7-2) 

d 2 u3'/dx 2 +d 2 u3'/dy 2 +d 2 uy/dz 2 =(w/E)(d 2 uy/dt 2 ), (7-3) 


a ] = [ cos{k x2 Li } 

_b\ [-sinfeAl 


where w'=w+ Aw is mass per unit length of the second region, 
with the TM attached. 

The vibrating amplitude solutions for the modified beam 
are now written u 1 '(x,y,z,t), u2’(x,y,z,t) and u3'(x,y,z,t) for the 50 
respective first, second and third regions. The boundary con- 

[ a 1 = T cos {k x2 Lz} 
b\ [sm{k x2 L2} 
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The vibration amplitude solutions are written as 


uY{x,y,z,t)=sm {k xi (-L x -x)] sin {k y y] sin {kyz) cos 
{co'o (t-t 0 )}. 

(12) 

u2'(x,y,z,t)={a sin {k x2 x}+b sin {k x2 *}} sin {k y y\ sin 
{Kz) cos {u>' 0 (t-t 0 )}. 

(13) 

u3'(x,y,z,t)=sm {k x3 (L-x)] sin (k ^ y } sin {k^) cos {co’ 0 
(‘-to)}, 

( 14 ) 

where k xl , k^ and k^ are wavenumber equivalents (possibly, 
all different) of the quantity k^ in Eq. (3), and a and b are 
parameters whose values are to be determined from the inter- 
face boundary conditions, Eqs. (8)-(ll). The frequency of 
variation oo' 0 for the first, second and third regions must be the 
same so that, in accord with the analog of Eq. (4) 

o# = (E/w){k 2 xl + k) + k\) = (EMV& + k 2 y + kjl = 

( 15 ) 

(E/WM&+AJ 

+ k 2 }. 

from which one infers that 


k xi =k x3 , 

( 16 ) 

k x2 -k x 1 2 = { ( w - w)/Z} 0) ' 0 2 . 

( 17 ) 

and the frequency oo'q can be detected by an appropriate 
mechanism, discussed in Appendix A. 

Equations (8) and 99) require that 

-sin {k x] [L x -L{)}=a cos {k x2 L{\-b sin {k x2 L { }, 

( 18 - 1 ) 

-k xl cos {k x3 (L x -L i ) } =k x2 { a sin {k^L^+b cos 

( 18 - 2 ) 


from which the parameters a and b are determined to be 


(1 /^sinfeLi } r -sin^ (L* - U )} (19) 

(1 / k X 2)cos{k x2 L>i } J [ —k x icos{—k x i (L^ - L { )} _ ' 


Equations (10) and (11) require that 

sin {k x3 (L x -L 2 )}=a cos {k x2 L 2 }+b sin {k x2 L 2 }, (20-1) 

-k x 3 cos {k x i {L x -L 2 ) } =k x2 {-a sin {k x2 L 2 }+b cos 

{K 2 L 2 }} (29-2) 

from which the parameters a and b are determined to be 


-(l/^sinfe^Jirsinfet^ ~Li)} (21) 

(1 / k x2 )cos{k x2 L 2 } J|_ -^ 3 cosfe(Z^ - L2)} _ ' 


ditions now include Eq. (2) and the interface relations (re- 60 Equating the parameters a in Eqs. (19) and (21) yields a first 
quired by continuity) dispersion relation for the wavenumbers k^ and k x2 in terms 

ul\x=-L l ,y,z,t)=u2'{x=-L h y,z,t), (8) of the lengths L : and L 2 . 

dul \x=-L h y,z, t)!dx =du2'(x=-L i ,y, z, t)/dx, (9) 

u2'{x=L 2 ,y,z, t)=u3 ' (x=L 2 ,yz, t ) , (10) ^ -k x 2 { cos {k^L^} sin {^(Z^-ZJj+cos {k x2 L 2 } sin 

{k x { (L x -L 2 ) } } =k x 1 { sin {k x2 L 2 } cos {k xl (L x -L 2 ) 

}+sin cos fc^-LT}}. 


du2'(x=L 2 ,y,z,t)ldx=duy{x=L 2 ,y,z,t)ldx. 


( 11 ) 


( 22 ) 
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Equating the parameters b in Eqs. (19) and (21) yields a 
second dispersion relation for the wavenumbers k^ (Ht^) 
and k x2 in terms of the lengths L l and L 2 . 

k x2 {sm {k x2 L { } sin {k^L^-L^-sm {k x2 L 2 } sin {k xl 
(L x -L 2 )}}=-k x i( cos {k x2 L 2 } cos {k xi (L-L 2 )}~ 
cos {k x2 L { } cos {k^-L,)}}. (23) 

For arbitrary values of L : and arbitrary values of L 2 , subject 
to the constraints -L X >-L 1 >L 2 >L X , and L 1 +L 2 =D (length of a 
single target molecule^constant), one can numerically deter- 
mine a minimum frequency difference 

Ao)(min)=min zl _ i2= ^l( 0 0 (ref)-(o' 0 1 (22) 

for presence of a single target molecule TM located at 
-L 1 ^x^L 2 . If an experimentally measured difference lco 0 
(ref)-co' 0 l is substantially less than Aoo(min), but substantially 
greater than zero, this indicates that at least one other mol- 
ecule, not the target molecule, ps present and is attached to the 
beam. Using this approach, one can distinguish between pres- 
ence of one or more target molecules and presence (only) of 
one or more other molecules, in some situations. 

More generally, when co' 0 (measured) differs substantially 
from oo 0 (ref) (measured or calculated), this indicates that at 
least one molecule (e.g., the target molecule TM) is present 
and is attached to the second region. The length D of the 
second region, on which the chemical substance CS is coated, 
is preferably no longer than the diameter (known) of the 
specified target molecule. Adopting this constraint will limit 
the molecule(s) that can become attached to the CS in the 
second regions so that, if a) , 0 (measured)=ci) , 0 (calculated, as in 
Eq. (15)), it is likely that the attached molecule is the TM. 
Further, only one TM can become attached to the single CS 
coating on the beam. 

EXAMPLE 

The chemical strip(s) 13 on the beam 11 in FIG. 1 contain 
a chemical substance CS, selected to serve as an quasi-anti- 
body, by binding to a target molecule (corresponding to an 
antigen) so that the antigen also becomes part of an extended 
beam, with a corresponding increased beam mass per unit 
area. A suitable example of the chemical substance CS is 
avidin, a tetrameric glycoprotein that occurs in egg white and 
in certain animal tissues, that has a combined mass of 67-68 
kiloDaltons (about 10 percent carbohydrate), and that binds 
tightly to, and inactivates, biotin. Biotin is a B vitamin and 
coenzyme that occurs, for example, in egg yolk. A variant of 
avidin, streptavidin, provides tighter and more specific bind- 
ing to bacteria and other biomolecules. Unlike avidin, strepta- 
vidin has no carbohydrate. Both avidin and streptavidin can 
be used to link molecules, such as radioisotopes and mono- 
clonal antibodies, together. Preferably, the avidin or strepta- 
vidin is deposited don the (bare) beam 11 as a monolayer, or 
in very few layers, and biotin molecules or biomolecules in 
the fluid 15 that pass by the beam may become strongly 
attached to the avidin or streptavidin. It may be preferable to 
use a gas, rather than a liquid, for some attachments and to use 
a relatively low transport rate of the fluid, to reduce the 
likelihood of detaching the biotin or detaching the avidin or 
streptavidin from the beam 11. Use of the avidin-biotin com- 
bination or the streptavidin-biotin combination is an example 
of the antibody-antigen combinations that can be used with 
the invention. 

An example of longitudinal-transverse beam vibration in 
an ultrasonic regime is presented by J. Tsujino and H. Yoshi- 
hara in “Vibration and Welding Characteristics of High Fre- 
quency Longitudinal -Transverse Vibration Ultrasonic Wire 
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Bonding Systems,” 1998 1.E.E.E. Ultrsonics Symposium, pp. 
735-738, where a steel rod, 7 mmx26.2 mm, with a free edge 
is found to have three longitudinal vibrating ranges, 351-368 
kHz, 433-472 kHz and 916-976 kHz, and to have transverse 
5 vibration ranges of 377-382 and 926-927 kHz. 

FIG. 3 is a flow chart of a method for practicing an embodi- 
ment of the invention. In step 3 1, one or more reference beam 
resonant frequencies co 0 (ref) are provided for a reference or 
unmodified beam that comprises a substrate, having one or 
10 more strips of a selected chemical substance CS that behaves 
as a quasi-antibody for a specified target molecule TM that 
behaves as a quasi-antigen. In step 32, a fluid, which may or 
may not contain the target molecule TM, is allowed to flow in 
a selected direction, transfers to a selected x-coordinate direc- 
15 tion, across the beam; and one or more target molecules, if 
present is allowed to attach to the strip(s) of chemical sub- 
stance CS on the beam, to provide a modified beam. In step 
33, one or more resonant frequencies o.y o , of the modified 
beam is estimated (e.g., measured experimentally), corre- 
20 sponding to the one or more resonant frequencies oo 0 (ref) of 
the reference beam. In step 34, the resonant frequencies, 
oo 0 (ref) and oo'o, are compared to determine if these two reso- 
nant frequencies are substantially the same. If the answer to 
the query in step 34 is “yes,” the system, in step 35, interprets 
25 this condition as indicating that the specified target molecule 
TM is likely not present, or is present in a very small concen- 
tration, in the fluid. If the answer to the query in step 34 is 
“no,” the system, in step 36, interprets this condition as indi- 
cating that the specified target molecule TM, or at least one 
30 other molecule that serves as a quasi-antigen, is present in the 
fluid and has attached itself to one or more of the strips of the 
chemical substance CS on the reference beam. 

Alternatively, a beam coating of length D may be posi- 
tioned on each of two opposed surfaces (e.g, top and bottom) 
35 of the beam, with each CS coating independently becoming 
attached to one of the target molecules TM. The analysis of 
this situation is similar to the situation discussed in the pre- 
ceding, except that the mass per unit length w' in Eq. (8) may 
have any of three values, wl (no TM attached), w2 (one TM 
40 attached or w3 (two TMs attached). 

Alternatively, a first beam coating of a first chemical sub- 
stance CS1 and length D1 may be positioned on a first beam 
surface, and a second beam coating of a second chemical 
substance CS1 (CS1^CS2) and length D2 (D1>D2) may be 
45 positioned on a second beam surface, where it is desired to 
estimate if a first target molecule TM1 (attached to CS1) or if 
a second target molecule TM2 (attached to CS2) is attached to 
the beam. It is assumed that ay o (calculated-TMl) differs 
substantially from co , 0 (calculated-TM2). Where oy o (mea- 
50 sured)=co' 0 (TMl -calculated) or where (jo f 0 (measured)=(o' 0 
(calculated-TM2), it is likely that the corresponding target 
molecule is present and attached to the beam. Where oo'o 
(measured)>max{co , 0 (calculated-TMl), ay o (calculated- 
TM2)}, it is likely that TM1 and TM2 are both present and are 
55 both attached to the beam. A separate calculation would be 
required to accurately estimate the frequency value oo'o (cal - 
culated-TMl and TM2 present), because of possible overlap 
of the x-coordinate regions for coating by CS1 and by CS2. 

60 APPENDIX A 

Estimation of Beam Resonant Frequency 

FIG. 4 schematically illustrates an embodiment for deter- 
65 mination of the wavenumber solutions {k x , k’ x }, or of the 
corresponding angular frequency solutions oo 0 , using diffrac- 
tion of electromagnetic signals. The beam is excited by Ion- 
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gitudinally vibrations at a selected angular frequency to, 
which may or may not coincide with a solution frequency to 0 , 
by a suitable excitation mechanism EM. A narrow beam, 
monochromatic electromagnetic signal 

ES& cos (A-l) 

having a selected excitation frequency co e , is directed at the 
beam at a selected incidence angle 0, relative to the longitu- 
dinal (x-coordinate) direction. 

The excited beam manifests a uniformly spaced sequence 10 
of undulations or ripples in the x-direction, illustrated in FIG. 

5, and a portion of the signal ES is reflected or diffracted from 
the undulating beam surface. One goal is to choose a combi- 
nation of excitation frequency and incidence angle (oo e ,0) for 
which the phase k e s-co e t for the excitation signal ES arriving 15 
at each peak of the undulation pattern differs by an integral 
multiple of 2 jt. Consider the two signal “rays” ES-1 and ES-2 
arriving at two adjacent undulation peaks, as in FIG. 4 (???). 
The distance between the two adjacent peaks is Ax=2jt/k x , 
and the distance from the location P2 to the location P3 on the 20 
undulation peak P3 is 

As 2 ^ = Ax sin 0. (A-2) 

and the time required for the signal ray to propagate this 25 

distance is 

As 2 ^ c= As 27 3kJ(a e , (A-3) 

where 

c=tojk e (A-4) 30 

is the velocity of electromagnetic propagation of the excita- 
tion signal ES in the ambient medium. Preferably, the elapsed 
time As 2 3 /c is an integer multiple K of a time increment 

35 

At(Pl,P3)=kJa 0 , (A-5) 

of a time required for a beam amplitude peak at position PI to 
move to the position P3; that is, 

As 2i kJ(H={2% sin Q)/(ck x )=KkJa 0 (K=l, 2,3,...), (A-6) 40 

where k x and oo 0 are determined as in the preceding analysis in 
the specification. 

Equation (A-6) determines the incidence angle 0 in terms 
of the known quantities k^, co 0 and c. When a match occurs, 45 
Eq. (A-6) allows an estimation of co 0 and/or k x for the vibrat- 
ing beam for a small integer multiple, such as K=l, 2 or 3. 
Experimentally, a resonant frequency co 0 (exp) of the beam 
can be estimated, using a thermocouple or similar heat mea- 
suring device and estimating an excitation angular frequency 50 
oo at which maximum beam heating occurs in a selected time 
interval, such as 5-30 sec. A resonant frequency can also be 
estimated by observing the beam vibration amplitude and 
estimating a frequency at which the amplitude is maximized. 

The estimated beam resonant frequency oo 0 (exp), after 55 
beam flow-by, is compared with the reference resonant fre- 
quency oo 0 (ref) for the “bare” beam in the absence of attach- 
ment of any target molecules or other quasi-antigens to the 
beam. Where the two resonant frequencies agree within a 
small difference, the system concludes that the target mol- 60 
ecule TM is not present in the fluid. Where the two resonant 
frequencies do not substantially agree, the frequency differ- 
ence can be used with the preceding analysis to estimate the 
mass per unit length w’ of a target molecule TM attached to 
the beam, from a measurements of resonant frequencies of a 65 
beam to which the target molecule TM is attached. To avoid 
ambiguity, it is preferable to provide an isolated beam coating 
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on a beam, as illustrated in FIGS. 6 A, 6B and 6C, with a 
length approximately equal to a single length D of the target 
molecule TM. 

APPENDIX B 

Estimation of Target Molecule Concentration 

Appendix A sets forth a method for experimental determi- 
nation of the wavenumber solutions or the corresponding 
angular frequency solutions for a single beam. Where two or 
more (preferably, three or more) beams, spaced apart in the 
fluid flow direction, are provided, as in FIGS 6A-6C, one can 
estimate concentration of the target molecule TM in the fluid. 

Consider an array of spaced apart beams 61 A-n (n=0, 1 , . . 

. , N) of length 2L X , as shown in FIG. 6A, where the beam 
center coordinate, x=x„, has a value, x„=x 0 +n xy (Ax=fixed 
positive value), where x 0 is a selected value. The distance, 
L x -L l5 of an interface between a beam portion that is coated 
with CS and an uncoated portion is preferably the same, but 
need not be the same, for each of the beams. Allow a fluid 
having a known concentration, r|=r| 1, of the target molecule 
to flow past the beams 61 A-n in the y-direction. Some, but 
probably not all, of the coated regions CR(n) will, as a result 
of the fluid flow-by of finite time duration, attach to and hold 
a target molecule TM; and as a result of such attachment the 
concentration of target molecules TM will decrease roughly 
exponentially as the downstream distance y increases. The 
target molecule concentration C(x,y,z) at a downstream dis- 
tance near a beam center can be approximated as 

C(x=x n , y,z)**C 0 exp{-an -Ay], (B- 1 ) 

where the parameter a (>0) is a measure of the affinity of CS 
for the TM and/or of an initial concentration of the TM, and an 
initial concentration C 0 is roughly constant («r| 1) throughout 
the fluid. The value of the parameter a will depend upon the 
concentration r|l and can be estimated by (1) performing a 
sequence of flow-by experiments for the same time duration 
on a sequence of the array of N beams shown in FIG. 6A; (2) 
allowing a first subset of the TM in the fluid to become 
attached to and held by a second subset of the coated regions; 
(3) measuring the cumulative number v(n)(=0, 1, . . . , n) of 
CS-coated regions that do not have a TM attached thereto; (4) 
repeating the steps (l)-(3) P times to obtain an average cumu- 
lative number v(n;mean) for the P experiments for the initial 
concentration r| 1; (5) forming an error estimate 


" 2 (B-2) 

£(ar) = y w n {v(n\ mean) - Coexp{-ar« ■ As}} ; 

H=0 


where w w are selected non-negative weight values and q is a 
selected positive value (preferably, q^l; for example, q=2); 
(6) estimating a=a(r|l) according to a minimization proce- 
dure, such as 

de(a)/da=0, (B-3) 

from which the value of a is estimated; and (7) repeating steps 
(l)-(6) for at least two other distinct initial concentration 
values, r|=r|2 and r|=r|3, for the fluid. 

The estimated value of a0(r|) will increase monotonically 
(approximately linearly over short intervals) with increasing 
TM concentration r|, except at very high concentrations, 
where the value of a0(r|) will likely saturate. It is preferable 
to perform steps (1 )-(6) for at least 3-10 concentration values, 



US 7,549,338 B1 


9 

spaced apart by multiplicative factors of 1 .5-3, to provide an 
interpolation curve, aO=IC (r|) versus r\, as illustrated in FIG. 

7. 

As an example of this minimization procedure, with q=2 
and 5 

w=exp(-aAAsr}, (B-4) 

where As is known by a not yet known, Eq. (B-3) can be 
re-expressed as a polynomial equation in u: 


a / \ /V 

d a = 2 ^ w n (n ■ As){(v(«; meari)w” - u 2n \ =0 ... 


When the (real) solutions u of Eq. (B-5) are determined or 
estimated, the corresponding values of a can be estimated 
from Eq. (B-4). Among a sequence of values of a, the smallest 
positive value is often preferred. Alternatively, the affinity 
parameter a can be estimated as 20 

a=l n{C(x=x i ,y,z)IC(?c=x n ,y,z)}l(x n -x ] ), (B-6) 

which an be used to initially estimate the TM concentration, 
using IC(a). 25 

To obtain an estimate of the (unknown) concentration r\ of 
the TM in the fluid of interest, this fluid is allowed or encour- 
aged to flow past the sequence of N beams, and each of the N 
beams is interrogated, for example, as in Appendix A, to 
determine which beams in the sequence have a TM attached 30 
thereto. The monotonic sequence v(n) is determined, as in 
step (3), and a value, a=aO, that minimizes the error sum €(a) 
in Eq. (B-2) is estimated. The interpolation curve aO=IC(r|) is 
then examined to estimate the corresponding TM concentra- 
tion value solution r|0 for which a0=IC(r|0). This value, 35 
r|=r|0, is an estimate of the TM concentration in the fluid of 
interest. 

The sequence or array of beams may also be arranged with 
longitudinal (x-) axes perpendicular to the fluid flow direc- 
tion, as in FIG. 6B, or with longitudinal (x-) axes transverse at 40 
an angle § relative to the fluid flow direction, as in FIG. 6C. 
What is claimed is: 

1. A method for sensing at least one or presence of, and 
concentration of, a target molecule that may be present in a 
fluid, the method comprising: 45 

providing a first thin beam, having a strip of a selected 
chemical substance of a selected length D, oriented sub- 
stantially parallel to a longitudinal axis of the beam, on 
at least one surface of the beam, where the selected 
chemical substance binds to a selected target molecule; 50 
providing a first resonant frequency col of vibration of the 
first beam and the at least one strip; 
allowing a selected test fluid, which may contain a target 
molecule in an unknown concentration, to flow past the 
first beam for a selected time duration At, to provide a 55 
modified first beam; 

providing an estimated second resonant frequency co2(est) 
of vibration of the modified first beam; 
providing a measured second resonant frequency 
oo2(meas) of vibration of the modified first beam, where 60 
oo2(est) is substantially different from col; 
when the frequency oo2(meas) is substantially equal to col, 
interpreting this condition as indicating that the target 
molecule is not present in the fluid; and 
when the frequency co2 (meas) is substantially equal to 65 
co2(est), interpreting this condition as indicating that the 
target molecule is present in the fluid. 
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2. The method of claim 1, wherein said process of provid- 
ing said first resonant frequency col comprises: 

providing a mathematical solution for a partial differential 
equation that expresses an amplitude u(r,t) of a longitu- 
dinal vibration of said first reference beam with at least 
one selected boundary condition imposed, where r rep- 
resents a vector of at least one location coordinate x for 
an arbitrary location on said first beam and t represents a 
time variable, where x is measured along said first beam 
longitudinal axis; and 

estimating said first resonant frequency col from a first 
dispersion relation that relates col to at least a first 
parameter k x that is determinable from the at least one 
boundary condition imposed on the solution u(r,t) for 
said first beam, where said target molecule is not present. 

3. The method of claim 2, wherein said process of provid- 
ing said estimated frequency oo2(est) comprises: 

providing a mathematical solution for a partial differential 
equation that expresses an amplitude u'(r,t) of a longitu- 
dinal vibration of said modified first beam with at least 
one selected boundary condition imposed, where r rep- 
resents a vector of at least one location coordinate x for 
an arbitrary location on said modified first beam and t 
represents a tine variable, where x is measured along 
said modified first beam longitudinal axis; and 

estimating said first resonant frequency co2(est) from a 
second dispersion relation that relates co2 to at least one 
wavenumber parameter k x that is determinable from the 
at least one boundary condition imposed on the solution 
u'(r,t) for said modified first beam, where said target 
molecule is assumed to be present and to be attached to 
said modified first beam. 

4 . The method of claim 3, further comprising: 

providing an estimate of a minimum frequency difference 

Aco(min)=min{ lco2(est)-col I }, where at least one of said 
target molecule is present on said modified first beam; 
and 

when an estimated frequency difference I co2— col I is sub- 
stantially less than the minimum frequency difference 
Aco(min), interpreting this condition as indicating that at 
least one molecule, other than said target molecule, is 
attached to said modified first beam. 

5 . The method of claim 1, further comprising selecting said 
chemical substance to comprise either avidin or streptavidin 
and choosing said target molecule to comprise biotin. 

6. The method of claim 1, comprising: 

(i) providing an ordered sequence of N-l supplemental 
beams, numbered m=2, . . . , N(Ni^3), each substantially 
identical to and spaced apart from and substantially par- 
allel to said first beam, numbered m=l , and each supple- 
mental beam having a strip of said selected chemical 
substance of length substantially equal to D on at least 
one surface, where supplemental beam .no.m is spaced 
apart downstream from said first beam by a distance of 
approximately (m-1 )As in a direction of flow of said test 
fluid; 

(ii) allowing said test fluid, having a known initial concen- 
tration, r|=r| 1 of said target molecule, to flow across said 
first beam and across the N-l supplemental beams to 
provide N modified first beams, numbered m-1 , . . . , N; 

(iii) determining, for the N modified beams, whether beam 
no. m' has said target molecule attached thereto; 

(iv) determining a number v(n) of beams, among the beams 
no. m-1 , . . . , n, that have said target molecule attached 
thereto; 

(v) estimating a value, a=a0(r|), for which an error sum 
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s{a) = ^ w„|k(«) -^ lexp-j-or-rt- As"}| 

n=l 

is minimized, where w„ is a selected non-negative value; 

(vi) repeating steps (ii)-(v) for said test fluid, having a 
known initial concentration, r|=r|2 («r|l) of said target 
molecule in said test fluid; 

(vii) repeating steps (ii)-(v) for said test fluid, having a 
known initial concentration, r|=r|3 (W|l, «r|2) of said 
target molecule in said test fluid; and 

(viii) estimating a relationship, aO=aO(r|), between the 15 
value aO and a concentration r\ of said target molecule in 
said fluid from the values aO(01), a0(r|2) and a0(r|3). 

7. The method of claim 6, further comprising: 

(ix) providing said first beam, numbered m=l, and said 
ordered sequences of N-l supplemental beams, num- 
bered m=2, . . . , N (Ni=3), each substantially identical to 
and spaced apart from and substantially parallel to said 
first beam, numbered m=l , and each supplemental beam 
having said strip of said selected chemical substance of 
length substantially equal to D on at least on surface, 
where supplemental beam .no. m is spaced apart down- 
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stream from said first beam by a distance of approxi- 
mately (m-l)As in a direction of flow of said test fluid; 

(x) allowing said selected test fluid, having said unknown 
concentration of said target molecule, to flow across said 

5 first beam and across said N-l supplemental beams to 
provide said N modified first beams, numbered m-1 , . . 
• , N; 

(xi) determining, for said N modified beams, whether said 
beam no. m' has said target molecule attached thereto; 

to (xii) determining a number v’(n) of beams, among said 
modified beams, no. m-1, . . . , n, that have said target 
molecule attached thereto; 

(xiii) estimating a concentration value r\' for said target 
molecule in said selected test fluid for which an error 
sum 


N 

s'(a0) w n \v\n) -^ / exp{-ar0(^ / )‘ w '^)l 

n = 1 

is minimized; and 

(xiv) interpreting this condition as indicating that said con- 
centration value for said target molecule in said selected 
25 test fluid is approximately r\\ 



